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P-i. 

We calculate the leading logarithmic QCD corrections to the decay B s — > 77 in the 
two Higgs doublet model with tree level flavor changing currents (model III) including O7 
r S^ . type long distance effects. Further, we analyse the dependencies of the branching ratio 

\ Br(B s — > 77) and the ratio of CP-even and CP-odd amplitude squares, R = \A + \ 2 /\A~\ 2 , 

on the charged Higgs mass m#± and the selected model III parameters £ U,D including the 
leading logarithmic QCD corrections. It is found that to look for charged Higgs effects, 
the measurement of the branching ratio Br(B s — ► 77) is promising. 
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1 Introduction 



Rare B meson decays constitute one of the most important classes of decays since they are 
induced by flavor changing neutral currents (FCNC) at loop level in the Standard Model (SM). 
Therefore they allow us to test the flavor structure of the SM and provide a comprehensive 
information about the fundamental parameters, such as Cabbibo-Kobayashi-Maskawa (CKM) 
matrix elements, leptonic decay constants, CP ratio, etc. These decays are also sensitive to the 
new physics beyond the SM, such as two Higgs Doublet model (2HDM), Minimal Supersym- 
metric extension of the SM (MSSM) [§, etc. 

Among the rare B meson decays, the exclusive B s — > 77 decay has received considerable 
interest in view of the planned experiments at the upcoming KEK and SLAC-B factories and 
existing hadronic accelerators, which may measure the branching ratios (Br) as low as 10~ 8 . 
Since the B s — > 77 decay has two photon system, it is possible to study the CP violating effects 
and it can be easily detected in the experiments by putting a cut for the energy of final 
photons ||. Further, this decay can give information about physics beyond the SM. 

The B s -> 77 decay has been studied in the framework of the SM (2HDM) (§) 
without QCD corrections. It is well-known, that the QCD corrections to the b — > S7 decay are 
considerably large (see H - [|HJ and references therein). Therefore, one can naturally expect 
that the situation is the same for the inclusive b — ► 577 decay. Recently, the analysis with the 
addition of the leading logarithmic (LLog) QCD corrections in the SM [fTT[| - |Tl3f] , 2HDM |L4 



and MSSM [15| has been done and the strong sensitivity to QCD corrections was obtained. 

The calculations shows that the Br of the decay B s — ► 77 is enhanced with the addition 
of the QCD corrections in the SM. The extension of the Higgs sector (model II in 2HDM) 
brings extra enhancement to the Br. However, the present theoretical results are at the order 
of (~ 1CT 6 ) are far from the experimental upper limit Hl6| 

Br(B s -> 77) < 1.48. 10" 4 . (1) 

In the present work, we study this decay in the framework of the two Higgs doublet model 
with three level flavor changing neutral currents (model III) including perturbative QCD cor- 
rections in the LLog approximation, by imposing a method based on heavy quark effective 



theory for the bound state of the B s meson [11 1. Further, we improve our calculations with the 



inclusion of long-distance effects through the transition B s — > ^7 — > 77, which we call 07-type, 
sec 



11 for details. In the analysis, we use the constraints coming from AF = 2 (F = K, D, B) 



mixing, the p parameter, the ratio R e b xp = Y(Z — > bb)/T(Z — > hadrons) (= 0.2178 ± 0.0011) 
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and the CLEO measurement of the decay B — > X s 7 ||17||) for the selected parameters of model 



III [0, and we get an extreme enhancement for the Br(B s — > 77). We also predict an upper 
limit for the model III parameter £^ using the present experimental restriction for the Br of 
the decay B s -> 77, eq. ([!]). 

The paper is organized as follows: In section 2, we give the LLog QCD corrected ampli- 
tude for the exclusive decay B s — > 77. We further calculate the CP-odd A~ and CP-even 
A + amplitudes in an approach based on heavy quark effective theory, taking the LLog QCD 
corrections into account. Section 3 is devoted to the analysis of the dependencies of the Br and 
the ratio |A + | 2 /|A~| 2 on the parameters /i (scale parameter), the Yukawa couplings £^ and 
our conclusions. In the appendix , we present the operator basis and the Wilson coefficients 
responsible for the inclusive h — > S77 decay in the model III. We also discuss the contributions 
of the neutral Higgs bosons to the coefficient C7 for the decay b — > 57 and the restrictions to 
the free parameters where one of the indices i or j belong to the first or second generation. 

2 Leading logarithmic improved short-distance contri- 
butions in the model III for the decay B s — > 77 

In this section we present the effective Hamiltonian to the exclusive B s — > 77 decay amplitude 
in the 2HDM with tree level neutral currents (model III). We start with the brief explanation 
of the model under consideration. The Yukawa interaction in the general case is 

Cy = vYjQiLfaUjR + v5QiL<PiD jR + ^Q iL ^U jR + gjQiL<hD jR + h.c. , (2) 

where L and R denote chiral projections L(R) = 1/2(1 =|= 75), 0j for i = 1, 2, are the two scalar 

u,d , t u r 
<ij and C, ; 

\ V2x 



doublets, rffj D and &' D are the matrices of the Yukawa couplings. The choice of <pi and fa 



* 1 = 7I 

with the vacuum expectation values 



v + H° ) + { t X ° 



>2 



J_( V2H+ 



(3) 



<^>=^[l) ;<0 2 >=o, (4) 

permits us to write the Flavor Changing (FC) part of the interaction as 

Cy,fc = &iLfoU jR + ^QiLfaDjR + h.c. , (5) 
where the couplings £ U,D for the FC charged interactions are 

£ci = ^neutral VqKM , 

ich — VCKM ^neutral , (6) 



and Neutral 1 * s defined by the expression 



Here, the charged couplings appear as a linear combinations of neutral couplings multiplied by 
Vckm matrix elements. 

Now, we would like to discuss the LLog QCD corrections to b — > 577 decay amplitude in 
the model III. As it is well known, the effective Hamiltonian method, obtained by integrating 
out the heavy degrees of freedom, is a powerfull one. In the present case, t quark, W ± , H ± , Hi, 
and H2 bosons are the heavy degrees of freedom. Here H ± denote charged, H\ and H2 denote 
neutral Higgs bosons. The LLog QCD corrections are done through matching the full theory 
with the effective low energy theory at the high scale /i = my/ and evaluating the Wilson 
coefficients from mw down to the lower scale /1 ~ 0(mb). Note that we choose the higher 
scale as n — m w since the evaluation from the scale \i = m H ± to fi = mw gives negligible 
contribution to the Wilson coefficients (~ 5%) since the charged Higgs boson is heavy due to 
the current theoretical restrictions based on the experimental measurement of B — > Xgj due 



to the CLEO collaboration |L7f, for example, m% > 340 GeV @, m% > 480 GeV (L§ 



The effective Hamiltonian relevant for our process is 



n eff = -^v tb v* £ dMOiin) , (8) 



where the Oi are operators given in eqs. (|23|), (j24[) and the Ci are Wilson coefficients renormal- 
ized at the scale /i. The coefficients Ci are calculated perturbatively. The explicit forms of the 
full-set operators and the corresponding Wilson coefficients are presented in Appendix. 

To obtain the decay amplitude for the B s — > 77 decay, we need to sandwich the effective 
Hamiltonian between the B s and two photon states, i.e. < B s \Ti. e ff\ / -f-f >. This matrix element 
can be written in terms of two Lorentz structures 0-0, JTl]]-[0: 



A(B S -> 77) = A + T^ v + iA-TpjFV , (9) 

where T^y = le^^JF ^. Here, the two different operator sets (eqs. (p3|) and (p4])) give 
contributions to the CP-even A + and CP-odd A~ parts. We denote the CP-even (CP-odd) 
amplitudes due to the first set as Af (Ai) and the second set as A^ (A2). In a HQET inspired 
approach these amplitudes are: 

A+ a em G F f Bs (I m% s {m e b Sf -m e Jf) ff 
1 " V2tt m\ s t V3A - (m B .-A a )(m? // + m e J f ) 7 KN 



1 In all next discussion we denote Neutral as £n' D ■ 
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9 ml + 



' ' m b J(m b ) + m s J(m s ))D(n) - m c J(m c )E(p,)' 



b t m s 



eff 



+ ^/ ef/ ^jj-[{m b A(m b ) +m s A(m s ))D(fi) + m c A(m c )E(fi)]] , (10) 

9(TO b JJ + m s J ) 



and 



+ a em G F f Bs x (I m Bs (m e b ff -m e / f ) , e// 

2 " v^tt m| s t \3A a (m B ,-As){ml ff + m e s f ^ 7 W 



^ = -^S^fnM (~ t-j4 ^:g-C?"to) , (11) 

V27T \3m B A s (m Bs -A s ) J 

where Q q = | for g = u, c and Q 5 = — | for q = d, s, b. Here, we have used the unitarity of the 
CKM-matrix J2i= u ,c,t V* s Vi b = and have neglected the contribution due to V* s V ub <C V t *V tb = 



\ t . The function g_ is defined as [1 1 



= m Bs ( m f + m f/) 2 + A s (m 2 Bs - (mf f + m e / f ) 2 ) . (12) 

The parameter A s enters eq. ( |T0| ) and ([11]) through the bound state kinematics (Tl[. Now, we 
would like explain the approach we follow to be clear how to get the amplitudes ( eqs. ( PH ) 
and flllD). The momentum of the b— quark inside the B s meson can be written as p = m b v + k, 
where k is the residual momentum and v is the 4- velocity defined by v = p B /ms s - Here p B is 
the 4-momentum of the meson B s . In the matrix elements we need to evaluate p.kj and p'.ki, 
i = 1,2, where ki, p' are the momenta of the outgoing photon and s -quark respectively. We 



average the residual momentum of the b quark as [20 



< >= -— !— (Aj + 3A 2 )f (U , 
2m b 



< k^k v >= y (g^ - v^vj) , (13) 

where Ai, A2 are matrix elements from heavy quark expansion. A2 (= 0.12 GeV) can be deter- 
mined from the B* s ^ — B^ mass splitting. 
Finally we have 

P-h = —{m Bs -A s ) , 
P-h = - T ^ ± A S , 

/ e//\2 2 2 o\ 

(m b JJ ) = p = m b - 3A 2 , 

( m ^//) 2 = p' 2 = (miff) 2 -rn 2 B + 2m Bs A s . (14) 



Here, we use the definitions Pb = P — p', — -yS v -k% = and the HQET relation |20| 



m Bs = m b + A s - ^-(\i + 3\ 2 ) , (15) 

In eq. flT4j), ml and m e /f are the effective masses of the quarks in the B s meson bound 
state and the parameter A s can be defined as A s = A + Am. A together with Ai can 
be extracted from data on semileptonic B ± , B decays and the measured mass difference 
Am = mg s — mg = 90MeV [J22| . Using eq. (JIB]), we correlate the parameters A s and m^ and 
A and Ai (Table (D). 

After this brief explanation of the HQET inspired approach we follow, we continue with the 
LLog QCD-corrected Wilson coefficients Ci...io(/f) Q] " Ell which enter the amplitudes in the 
combinations 

cm = cm = {cm - cm)n c + cm - cm , 
cm = (cm + cm - cm)n c + cm + CM - CM , 
CM = cm = (cm + cm){n c + i) - n c cm - cm , 
D(fi) = cm + cmn c , 

E(fi) = C x M + CMNc (16) 

where N c is the number of colours (N c = 3 for QCD). 

The effective coefficients Cj^(^l) and Cy •'(/i) are defined in the NDR scheme, which we 



use here, as 18 



C e 7 ff (v) = Cr DM (n) + Q d (Cl HDM (n) + N c Cl HDM ^)) , 

+ Q u (^ CIq DM (ii) + N c ^ Cr DM {^)) , 
m b m b 

C' 7 eff (») = C' 7 WDM ^) + Q d {C^ HDM ^) + N c C'i HDM ^)) 

+ Q u {^ C? HDM (n) + N C ^ Cr DM (»)) ■ (17) 
m b m b 

The functions I(m q ), J(m q ) and A(m g ) come from the irreducible diagrams with an internal q 
type quark propagating and are defined as 

2 

771 

I(m q ) = 1 + A(m 9 ) , 

m B s 

m n — 4m? 



A(m n 



2 arctan( 



Am? — m 



— 7r for 



m 



B. 



4m^ 



Finally the CP-even A and CP-odd A + amplitudes can be written as 



< 1. 



A + = A+ + A+ 
A' = A[ + A^ 



(19) 



where the amplitudes Af and Af are given in eqs. ( fLOD and (pi] ) respectively. 
In our numerical analysis we used the input values given in Table ([!]). 



Parameter 


Value 


m c 


1.4 (GeV) 


m b 


4.8 (GeV) 


a- 1 

"em 


129 


At 


0.04 




4.09 ■ 10~ 13 (GeV) 


IBs 


0.2 (GeV) 


m Bs 


5.369 (GeV) 


m t 


175 (GeV) 


m w 


80.26 (GeV) 


m z 


91.19 (GeV) 


A (5) 
IV QCD 


0.214 (GeV) 


a s (m z ) 


0.117 


A 2 


0.12 (GeV 2 ) 


Ai 


-0.29 (GeV 2 ) 


A s 


590 (MeV) 


A 


500 (MeV) 



Table 1: Values of the input parameters used in the numerical calculations unless otherwise 
specified. 
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3 Discussion 



Before we present our numerical results, we would like to discuss briefly the free parameters of 
the model under consideration, i.e. model III. This model induces many free parameters, such 
as ^fj D where i,j are flavor indices. To obtain qualitative results, we should restrict them using 
the experimental measurements. The explicit expressions of C7 and Cf° (eq. (g§) show that 
the the neutral Higgs bosons can give a large contribution to the coefficient C-j which can be 



in contradiction with the CLEO data 17 



Br{B -> X s7 ) = (2.32 ± 0.07 ± 0.35) 10" 4 . (20) 

Such potentially dangerous terms are removed with the assumption that the couplings ^ is (i = 
d, s, b) and £jy db are negligible to be able to reach the conditions £jy bb ^ is « 1 and ^ db £® ds « 



1. Further, we use the constraints [18| coming from the ratio R e b v = Y{Z — ► bb)/T(Z — > 
hadrons), namely £®bb > ^p 2 ^ the restrictions due to the AF = 2 mixing, the p parameter 
J23|, and the measurement by CLEO collaboration. The analysis of the mentioned processes 
and the discussion given above, leads to choose ^tc << ^ tt , £§ bb and £® ib ~ , ~ 0> where 
the indices i,j denote d and s quarks . 

After these preliminary remarks, let us start with our numerical analysis. In this section we 
study the dependencies of the Br and the ratio R = \A + \ 2 /\A~\ 2 on the selected parameters 
of the model III (£n U , ^bb an d w#±) and the QCD scale \i. In figs. [l| and || ( || and |J) we 
plot the Br of the decay B s — > 77 with respect to the charged Higgs mass for the fixed 
value of £,Nbb = 60 m b (£,§ bb = 90 m b ) at three different /i scales (/i = 2.5 , 5 , m w )GeV . Fig. |l] 
(H) represents the case where the ratio \r t b\ = |*zr^ <<: ^ anc ^ snows that the Br obtained 
in the model III almost coincides with the one calculated in the SM. The suppression of the 
contribution coming from the charged Higgs boson can be seen from eqs. ( P5| ) and |IT] with the 
choice \r t b\ « 1- However in fig. ||] @, the Br is presented for the case r#, >> 1. It is seen 
that there is an extreme enhancement of the Br, especially for the small values of m#±. Note 



that this is similar to the result coming from the choice of tan/3 < 1 in the model II 14]. In 



figs. |5| (H) we present dependence of the Br at the fixed value of m H ± = 500 GeV for 
\r tb \ « 1 {r tb » 1). In the region \r tb \ « 1, the Br is nonsensitive to the coupling ^ bb 
and almost coincides with the SM value, however for r t b » 1, it increases with the increasing 
i^bb- Note that the Br is sensitive to the scale /1. For \r tb \ « 1, it increases with decreasing 
\x (figs. H and |] ) similar to the model II Ul4| . For r tb >> 1, decreasing the scale /i causes 
the Br to decrease, (figs. |], § an d @). 
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At this stage we would like to estimate the upper limit of £j} bb for r tb » 1 by using the 
present experimental result Br(B s — > 77) < 1.48 • 10~ 4 . By choosing the lower limit for the 
mass m H ± ~ 480 GeV JTJJ at the scale /1 = 2.5 GeV, we get ^ bb < 90 m b . It is interesting to 
note that the increasing value of m#± makes the restriction region for £jy bb smaller. 

Since the two photon system can be in a CP-even and CP-odd state, B s — > 77 decay allows 
us to study CP violating effects. In the rest frame of the B s meson, the CP = — 1 amplitude 
A~ is proportional to the perpendicular spin polarization e[ x e^, and the CP = 1 amplitude 
A + is proportional to the parallel spin polarization eLe^. The ratio R is informative to search 
for CP violating effects in B s — > 77 decays and it has been studied before in the literature in 
the framework of the 2HDM without || and with [[TJ| QCD corrections. 

Now, we present the dependence of the ratio R = \A + \ 2 /\A~\ 2 on the selected parameters of 
model III, displayed in a series of figures ( [7| - [□]). In figs. |7| and ^| we plot the dependence of R 
on rriH± for fixed £® bb = 60 m b and three different \i scales, (2.5 , 5 , mw) GeV. The R ratio is 
almost nonsensitive to m H ± for \r tb \ << 1 (fig. |7]). However, it is enhanced with the increasing 
value of m H ± {m H ± < 1000 GeV) for r tb » 1 (fig. |8|). Decreasing the scale /i weakens the 
dependence of the ratio R on m#± and the contribution of the charged Higgs bosons to the R 



ratio becomes small similar to the model II [14|. Further, it becomes less dependent to 
with increasing ^ bb (fig g, g). 

Fig. |H] and |ll] show the dependence of R on ^ bb for fixed vtlh± = 500 GeV. Like the 
previous case, R dependence to the coupling ^ bb is extremely weak for \r tb \ « 1. (fig. |T0|) . 
For r tb » 1 the ratio R increases with decreasing ^ bb (fig [11]). Note that, this ratio can 



exceed one unlike the SM case. The same situation appears in model II also [Tj] , 

The ratio R is quite sensitive to QCD corrections and it is enhanced with decreasing scale \x 
for the SM and \r tb \ « 1 (figs. [7], |TDD . However, this ratio decreases for r tb » 1 (figs. ||, |9|, |TT| ). 
We observe, that the smaller the value of m H ± (the larger the value of ^ bb ), the less dependent 
is the ratio on /1. This strong fi dependence makes the analysis of the model III parameters 
m#± and £ U ' D for the given experimental value of the ratio R to be difficult, especially for the 
case \r tb \ << 1. However, we believe, that the strong \i dependence will be reduced with the 
addition the of next to leading order (NLO) calculation. This requires a computation of finite 
parts of many two loop diagrams, and divergent part of three-loop diagrams. This lies beyond 
the scope of the present work and it has to be done to reduce the uncertainity coming from the 
scale \i. Fortunately, the choice of fi = m b /2 in the LLog approximation reproduces effectively 
the NLO result for the b — > S7 decay and one can suggest that it may also work for the b — > S77 
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decay. In any case, the analysis on the model III parameters will be more reliable after NLO 
calculations are done. 

We complete this section by taking the 7 type long distance effects (LD 0r ) for both the 
Br and the ratio R into account. The LDq? contribution to the CP-odd A~ and CP-even A + 
amplitudes has been calculated with the help of the Vector Meson Dominance model (VMD) 
|TT| and it was shown, that the influence on the amplitudes was destructive. With the addition 



of the LD 0t effects, the amplitudes entering the Br and R ratio are now given as 

A- = A- SD + A- LDo7 , (21) 

where A^ D are the short distance amplitudes we took into account in the previous sections 
(eq.flilf)). The LDq 7 amplitudes A^ Do ^ are defined as |JTT 



4+ _ K a emG F - m b {m Bs - m ) ff 

ai D07 = v^^^mumt^-c^^) , (22) 

where /^(0) = 0.18 GeV is the decay constant of <j) meson at zero momentum, Fi(0) is the 
extrapolated B s — > form factor (for details see [fLl]| ). Note that we neglect the contribution 
of operator 0' 7 in eq.( |22|) since the coefficient C' 7 ' e (n) is negligible compared to the coefficient 
C?^(ji) (see the discussion given in appendix). 

In figs. ( [H - O ) we present the m H ± and bb dependencies of the Br and the ratio R 
with the addition of LDo 7 effects. Here we use ^i(O) = 0.16 The Br decreases with the 
addition of LDq 7 effects, since the effect is destructive. The /1 scale uncertanity of the Br is 
smaller compared to the case where LD effect is not included. 

It can be shown, that the value of the ratio R also decreases with the addition of LDo 7 
effects for the SM case. (figs. [15]- |T7|). However, while m#± is increasing or £jy 66 is decreasing, 
the effect of the LDo 7 contribution on the ratio R is increasing for r tb » 1 (figs. [TB|, |X7|) . 

There are still non-perturbative effects which can come from the formation of cc bound 
states. However these states are far off-shell and do not give significant contribution to the 



decay rate [JT3] . For example the chain process B s (f)ip — > 07 — > 77 is estimated and it is 



found to be at most 1% of the branching ratio Br(B s ^^)sd+ld 

Besides the strong [i dependence there is another uncertainity coming from the choice of 
bound state parameters m b and A s . It follows that the larger m b (smaller A s ), the larger 
Br and R ratio. Here the enhacement of the Br is caused by the j- dependence in amplitudes. 

9 



In conclusion, we analyse the selected model III parameters ( £iv&6> Cnu> m H ± ) an d QCD 
scale /i dependencies of the Br and R ratio for the decay B s — > 77. We predicted the upper 
bound for £j} bb , ^Nbb — 90 m b in the case of r tb » 1, using the constraints for £,® bb , £jy tt |18 
and the experimental upper limit for the Br of the decay underconsideration. We obtain that 
the strong enhancement of the Br is possible in the framework of the model III. 

A Appendix 

The operator basis and the Wilson coefficients for the 
decay b —> 577 in the model III 



The operator basis is the same as the one used for the b — > 37 decay in the model III |18| and 
SU(2) L x SU(2) R x U{1) extensions of the SM |5|: 

Ox = {s La l^c Lf3 )(c L ^b La ), 

2 = {sLal^C La )(c L ^b L p), 

3 = (s La -ff,b La ) Y (VLpYgLp), 

q=u,d,s,c,b 

A = (s La 7 M 6 L/3 ) Y {QLffQLa), 
q=u,d,s,c,b 

q=u,d,s,c,b 

6 = {sLal^bip) Y {QR^YQRoc), 
q=u,d,s,c,b 

e 



°7 = T^? s a (J ^( m bR + r n s L)b a J rfll/ 
9 = (s La 'j^c L f 3 )(c R p'j tl b Ra ), 



O10 = {s La l^c La )(c R p^b R p), (23) 
and the second operator set 0[ — 0[ Q which are flipped chirality partners of 0\ — Oi : 

0[ = (s R a r r l MCR/3)(cR l3 J fl bBa), 

0' 2 = {s Ra -i^c Ra ){c R ^b R p), 

q=u,d,s,c,b 

O4 = {s R al^b R p) Yj (QWfQRa), 
q=u,d,s,c,b 

0' 5 = {sRalpbRa) Y (QLpfqLp), 
q=u,d,s,c,b 
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0' 6 = (sRa^bnp) (QLprfqLa), 
q=u,d,s,c,b 

6 

°* = -^2^T^a^{m b L + m s R)b^\ 
0' g = (s Ra j^c R p)(c L ^b La ) , 

O' l0 = (s Ra 'y ll c Ra )(c Ll rfb L i3) , (24) 

where a and (3 are £77(3) colour indices and T^ v and Q^ u are the field strength tensors of the 
electromagnetic and strong interactions, respectively. 

In the calculations, we take only the charged Higgs contributions into account and neglect 
the effects of neutral Higgs bosons. At this stage we would like to give the reasons by using the 
restrictions to the effective Wilson coefficient Cy coming from the CLEO data in the process 
B — > K*j. (see section 2 and |L8| also) 

The neutral bosons Hq, Hi and H% are defined in terms of the mass eigenstates Hq ,ho and 
A as 

H Q = (H Q cosa — h sina) + v , 
Hi = (h cosa + H sina) , 

H 2 = A , (25) 

where a is the mixing angle and v is proportional to the vacuum expectation value of the 
doublet 0i (eq. (Q)). Here we assume that the massess of neutral Higgs bosons h and A 
are heavy compared to the b-quark mass. The neutral Higgs scalar ho and pseduscalar Aq 
give contribution only to C7 for our process. With the choice of a = 0, and can be 
calculated at m-w level as 



c^( mw ) = {v tb v: s y l £ il bi I 



D tD 

NAs 



Cf°(m w ) = {V tb V t :r l J2 ®MtS,i.a^—, (26) 

where and Qi are the masses and charges of the down quarks (i = d, s, b) respectively. Here 
we used the redefinition 



= f-^ ^ D • (2T) 

Eq. shows that neutral Higgs bosons can give a large contribution to C7, which does not 
respect the CLEO data JT^- Here, we make an assumption that the couplings ^§ is (i = d,s,b) 
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and Civ* are negligible to be able to reach the conditions £tv,&6 £iv,is 1 an d ^N,db^N,ds <<: 1- 
These choices permit us to neglect the neutral Higgs effects. 

Denoting the Wilson coefficients for the SM with Cf M (m\y) and the additional charged 
Higgs contribution with C^(mw), we have the initial values for the first set of operators 
(eq.( p3]) ) ( fll"8| and references within) 



C l,3,...6,9,lo( m W r ) - , 

C s 2 M {m w ) = 1, 

„a M . , 3x 3 -2x 2 1 -8a; 3 - 5x 2 + 7x 



4(x - l) 4 24(x - l) 3 

^"w - - "" + st i)t ■ 

Cl,..6,9,lo( m w) = , 

c?W) = A(^,« + ^,*c^)(S.« + S. te ^)J f i(y), 

T^t ^is V tb 

™t m b V ts V t b 

m t Vts V tb 

+ ^c^i) + £5U^)G a (y) , (28) 



and for the second set of operators eq. 



C'{ SM io(m w ) = 

C l H ...6,9,w( m w) 

C' 7 H (m W ) = —2 (^N,bs7J^ + ^N,ss) i^N,bb + ^N,sb7T-) F l(y) ^ 
m + vts v t b 



m t m b 'Vts v tb 

C' 8 H (m w ) = —2 {^N,bs^ + (£at,6& + ^N,sb7r)^{y) > 

+ — (^ >ta ^ + J + S, J^Ci/) , (29) 
where x = mf/m^ and y = mf/m 2 H ±. The functions F±(y), F 2 (y), Gi(y) and G 2 (y) are given 



as 



72(y-l) 3 12(y-l) 4 y 
y(5y - 3) , y(-3y + 2) ; 
Fa(v) = T%^l7 + 6(,-l) 3 ,ny ' 
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G >(w) = 4(^iy + 2(^Ty lny - (30) 

Note that we neglect the contributions of the internal u and c quarks compared to one due to 
the internal t quark. 

For the initial values of the Wilson coefficients in the model III (eqs. (^)and (p^)), we have 

Cl HDM {m w ) = 1, 

C 2 7 HDM (m w ) = C? M {m w )+C?{m w ) , 
Cl HDM {m w ) = Cl M {m w ) + C» {m w ) , 

r i2HDM I \ _ n 
U l,2,3,...6,9,W\ m W) - U , 

C' 7 2HDM (m w ) = C' 7 SM (m w ) + C 7 H (m w ) , 

Cf™(m w ) = Cf M (m w ) + Cf(m H ,) . (31) 

At this stage it is possible to obtain the result for model II, in the approximation ^ ~ 
and ~ 0, by making the following replacements in the Wilson coefficients: 

S*S = (32) 

and taking zero for the coefficients of the flipped operator set, i.e C[ — > 0. 

The evaluation of the Wilson coefficients are done by using the initial values Cf HDM 
(C'f HDM ) and their contributions at any lower scale can be calculated as in the SM case WB\ . 
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10*7 Br 



Figure 1: Br as a function of the mass m#± for fixed ^N,bb = 60 mj, at the region \r t b\ « 1. 
Here solid lines correspond to the scale \x = 2.5 GeV, dashed lines to /i — 5 GeV and small 
dashed lines to /x = GeV. The lines corresponding to the SM coincides with the lines we 
present here. 




Figure 2: The same as Fig 1, but at the region r t b » 1. Dotted dashed lines correspond to 
the SM. Note that dotted dashed line almost coincides with the ran axis 
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10*7 Br 



Figure 3: The same as Fig 1, but for fixed ^Nbb = 90 m& value 




Figure 4: The same as Fig 2, but for fixed £,jy bb = 90 m b value 



Figure 5: Br as a function of the coupling t;® bb for fixed m H ± = 500 GeV at the region 
\r tb \ « 1. Here solid lines correspond to the scale n = 2.5 GeV, dashed lines to /i — 5 GeV 
and small dashed lines to jj, = mw GeV. The lines corresponding to the SM coincides with the 
lines we present here. 
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300 400 500 600 



Figure 6: The same as Fig 5 but at the region r tb » 1. 



I A+l"2/ | A- I -2 



Figure 7: R ratio as a function of the mass 7n#± for fixed £jvw> = 60 at the region |r#,| << 1. 
Here solid lines correspond to the scale \i = 2.5 GeV, dashed lines to fx = 5 (7eV and small 
dashed lines to fx — m w GeV. The lines corresponding to the SM coincides with the lines we 
present here. 



I A+l"2/ | A- I -2 



1000 H 



Figure 8: The same as Fig 7, but at the region r tb » 1. Here solid curves (lines) correspond 
to the scale fx = 2.5 GeV, dashed curves (lines) to fx — 5 GeV and small dashed curves (lines) 
to fx = m w GeV for model III (SM). 
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I A+l "2/ | A- I -2 



Figure 9: The same as Fig 8 but for fixed £^ 66 = 100 m b 



I A+ | "2/ | A- | A 2 



Figure 10: R ratio as a function of the coupling C,^ bb for fixed m#± = 500 GeV at the region 
\r t b\ « 1- Here solid lines correspond to the scale \i = 2.5 GeV, dashed lines to /i — 5 GeV^ 
and small dashed lines to /x = Gel 7 . The lines corresponding to the SM coincides with the 
lines we present here. 



A+|"2/ | A- I A 2 



200 300 400 500 



700 800 



Figure 11: The same as Fig 7, but at the region r tb » 1. Here solid curves (lines) correspond 
to the scale /i = 2.5 GeV, dashed lines to /j, — 5 GeV and small dashed lines to /j, — m w GeV 
for model III (SM). 
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Figure 12: The same as Fig. 2 but including LD effects. The lines corresponding to the SM 
coincides almost with the m H ± axis. 
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Figure 13: The same as Fig. 5 but including LD effects. 
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Figure 14: The same as Fig. 6 but including LD effects. 
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A+|"2/ | A- I A 2 



Figure 15: The same as Fig. 7 but including LD effects. 



A+|"2/ | A- I -2 
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Figure 16: The same as Fig. 8 but including LD effects. 
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Figure 17: The same as Fig 11 , but including LD effects. 



21 



